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Abstract. - We analyze the dissipation of the vibrations of nano-mechanical devices. We show 
that the coupling between flexural modes and two-level systems leads to sub-ohmic dissipation. 
The inverse quality factor of the flexural modes of low frequencies depends on temperature as 
providing a quantitative description of the experimental data. 



Introduction. — Nano-electro- mechanical devices [1, 
2] (NEMS) are systems with great potential for applied 
' physics and engineering because of their extreme sensitiv- 
ity, as probes, to their environment [3-9]. Furthermore, 
because of their small size and large surface to volume 

■ ratio, these systems are in the crossover region between 
' classical and quantum behavior, and hence of great the- 
. oretical interest. Thus, the study of the sources of noise 
' and dissipation in these systems has attracted a great deal 
I of attention [10-18]. One of the common realizations of 
' nano-mechanical resonators is a rigid beam of nanoscopic 
. dimensions which vibrates at GHz frequencies [19,20]. The 
' damping of these oscillations has been a subject of intense 
I investigation [10, 13, 15, 16], as it sets a limit to their pos- 
sible applications. 

The damping of a low frequency oscillation in a NEMS 
I comes from the coupling of this mode to other low en- 
i ergy degrees of freedom. Experiments suggest that sur- 
' faces of nano-mechanical resonators resemble amorphous 

■ bulk systems [21], with a high density of defects, playing 
' a major role as a source of dissipation [22]. In amorphous 

solids, disorder and impurities lead to the existence of an- 
harmonic excitations, which can be modeled as a degree 
of freedom tunneling between two potential wells [23]. A 
simpler two-level system (TLS) description arises at low 
temperatures, when only the two lowest eigenstates have 
to be considered. Some properties of the distribution of 
TLSs in terms of their parameters (bias, Ag, and tunnel- 
ing rate, Ag) can be inferred from experiments [24], and 
they are considered to be the main source of damping of 
acoustical modes in disordered insulating solids [24-27]. 



In this work, we study mostly a rigid beam geometry, 
sketched in fig. [1], and analyze the dissipation processes 
for low energy flexural modes due to the presence of effec- 
tive TLSs at its surface. The generalization to dissipation 
of torsional modes is straightforward and will be given 
elsewhere. Once a given mode is externally excited, the 
TLSs living at the surface of the beam, which are coupled 
to this mode, will absorb part of its energy. But as the 
TLSs are also coupled to the rest of vibrational modes of 
the beam, they will release most of this energy to them. 
Thus the TLSs give rise to an indirect coupling between 
the externally excited mode and the rest of modes. For 
the experimentally relevant case of low amplitudes of vi- 
bration this coupling prevails over the usual anharmonic 
coupling. 

This energy flow process will be described in two stages: 

i) From the point of view of a given TLS, its coupling 
to the vibrational modes of the beam, which can be seen 
as an external bath, alters its dynamics and enables it 
to absorb and emit energy in a broad range of frequen- 
cies. In particular, the presence of flexural modes leads 
to the possibility of qualitative changes in the dynamics 
of the TLSs, as the former constitute a sub-ohmic envi- 
ronment [28, 29] for the TLSs. This is a consequence of 
the quadratic dispersion relation characteristic of flexural 
modes, which results in an enhancement of the density of 
low frequency modes. 

ii) Coming back to the externally excited vibrational 
mode whose damping we want to compute, the TLSs, 
dressed by all the vibrational modes of the structure, con- 
stitute the dissipative environment for the mode. This 
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Fig. 1: Sketch of the NEMS: a) Doubly clamped beam, sus- 
pended by both ends; b) Cantilever, with one free end. The de- 
vice is characterized by its width (to), thickness (€), and length, 
(L), where w ~ £ <C i. 



involves TLSs that are nearly resonant with the mode un- 
der consideration but also off resonance TLSs since the 
strong phonon/TLS coupling provides each TLS spectral 
function with tails far from resonances. In many NEMS 
experiments, the flexural modes studied are highly excited, 
either because of an external driving mechanism, or be- 
cause the temperature is much higher than the frequency 
of the mode. As a given TLS can absorb and emit over a 
broad range of energies, due to the incoherent tails in its 
spectrum (see below), these processes allow for the trans- 
fer of energy from a highly excited low frequency flexural 
mode to other modes with higher frequencies. 

Dissipative mechanisms other than TLSs, which have 
been extensively studied elsewhere [30-32], are not con- 
sidered here. 

TLSs coupled to low dimensional vibrations. 

The TLSs are mainly coupled to the strain induced by 
phonons. It is assumed that the main effect of the strain 
is to modify the energy splitting between the TLSs energy 
levels [33]. The Hamiltonian of a given TLS is character- 
ized by a bias, Aq, and the tunneling rate, Ag (we use 
units such that ?i = 1 = fc^, and omit the part of the 
Hamiltonian describing free vibrations): 



AqCTx 



(1) 



where diUj is a component of the deformation gradi- 
ent matrix, and F is an arbitrary function. Changing 
basis to the energy eigenstates of the TLS, eq.(l) be- 
comes H ^ Apg, -I- [( A|^/Ao)g:, -f (A^/Ao)f7;.]F(a^Uj). 
Ao = ^(Ag)2 ^ (Ag)2 is the splitting of the TLS. Dis- 
sipation is dominated by slightly biased TLSs for whom 
Ag <C Aq, so the last term can be ignored. A further 
expansion of F to lowest order in the displacement, to- 
gether with a 7r/4 rotation of the eigenbasis, leads to: 
i? = AqcTx + 7(Ag/Ao)crz9iUj , where 7 is the couphng 
constant (with dimensions of energy) . 

The main interaction between phonons and TLSs is due 
to the coupling to the operator cr^ of each TLS. Hence, the 
absorption properties of each TLS can be characterized by 



the spectral function: 

A(c^) =^|(0|a,|n)|'(5(w-w„+^o) (2) 

n 

where |n) is an excited state of the total system TLS plus 
vibrations. The linearization of the coupling implies that 
the interaction between a given TLS and the vibrations 

can be written as Rint = <^zJ2k (bl + bkj ■ 

We also define a spectral function that determines 
the damping induced by phonons on the TLS: J{uj) = 
12k I'^fcl ^ ^k) where ujk is the energy of mode k. 
The vibrational modes of a beam with fixed ends have a 
discrete spectrum, but we will approximate them by a con- 
tinuous distribution. This approximation will hold as long 
as many vibrational modes become thermally populated, 
kT ^ hujfund, where uifund is the frequency of the lowest 
mode. The condition is fulfilled in current experimental 
setups. 

Acoustic modes of a nanoscopic beam. Using contin- 
uum elasticity theory [34], a one-dimensional (ID) rod has 
compression and twisting modes, with a linear relation be- 
tween frequency and momentum, and bending, or flexu- 
ral, modes, where the frequency depends quadratically on 
momentum. We will consider a rod of length L, width w 
and thickness t, see Fig.[l]. We describe next the spectral 
function which describes how the modes of the rod absorb 
energy in different energy ranges. 

The compression and twisting modes lead to an ohmic 
spectral function for uj ^ 2ttc/R {R being a typical 
transversal dimension of the rod and c the sound veloc- 
ity), when the rod is effectively ID. In terms of the Young 
modulus of the material, -E, and the mass density, p, we 
get: Jcomp(^) = ad^l, where, 

= {^A-/Aof{2n^ptw)-\E/p)-y^ . (3) 

The twisting modes are defined by the torsional rigidity, 
C = ixt^w/i {fi is a Lande coefficient), and I = J dSx^ = 
t^w/\2 (where S is the cross-section). The corresponding 
spectral function is given by: Jtorsion(w) = atlw], where 

at = CijA-JAofi8n^ptwpI)-\pI/Cf^' . (4) 

The analysis of the flexural (bending) modes differs 
substantially from the other ones, because they corre- 
spond to two fields ^j{z,uj) {j = x,y) that satisfy [34]: 
Eljd^^j — ptw uj^ , where, for the system considered 
here, Ij — t^w/12. The normal modes have a quadratic 
dispersion ujj{k) — Elj / {ptw) k^. Their corresponding 
spectral function is sub-ohmic [28], Jficx{(^) — ceb^/uJ^^/uj, 
with. 



■=0.3- 



(l + i^)(l-2l/) //9\l/4 



S(3- 5i/) 



(I) 



(5) 



where v is Poisson's ratio and Wco — \/EIy/ {ptw){2T:/t)^ 
is the high energy cut-off of the bending modes. Collecting 
the previous results, we find the spectral function J{lu) 
plotted in Fig. [2]. 
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Fig. 2: Sketch of the contributions to the spectral function 
which determines the dynamics of the TLS. 

Dynamics of the TLSs. — We are interested in 
the TLSs that affect most the low-energy flexural vibra- 
tions. Hence, we focus on the TLSs whose tunneling 
amplitudes lie in the region where the damping is sub- 
ohmic. We define as the tunneling amplitude of a given 
TLS, including the renormalization due to the high en- 
ergy acoustic modes. We assume, as in the case of glasses, 
that the distribution of these TLSs is given by [24, 25] 
£((Ar, A^) = P/Ar- The sub-ohmic coupling, eq.(5), leads 
to a renormalization of A^: 

Aren = Arexp{-ab^/uJ^ J du) J {u) / Oj"^} (6) 

This equation has no solutions other than A^en = if 
A^ <C a^Wco, so that the tunneling amplitude of the low 
energy TLSs is strongly suppressed [28,29,35,36]. The 
remaining TLSs experience a shift and a broadening of 
the spectral function function A{ui), defined in eq.(2). In 
addition, acquires a low energy tail, which, at zero 

temperature, is [37]: 



A{w) CX ab UJ < Aren (7) 

ren 

There is also a high energy part, cx 
abV^coA^g^w""^/^, for Lo ^ Aren- The main fea- 
tures of A{u!) are shown in fig. [3]. Finally, we obtain the 
width of the resonant peak, T{Aren), using Fermi's golden 
rule, r(Aren) = 16af,^Wco V Aren- This description is 
valid for wavelengths such that, 1/L <C fc <C l/max(ti;, t). 

The value of A 2 is not renormalized by the phonons, so 
that the TLS cannot exchange energy with the environ- 
ment at frequencies lower that A^. In the following, we 
will consider only TLSs with Az<Aren- 

The total absorption rate by the ensemble of dressed 
TLSs present in the beam, Atoti'-^), is obtained summing 
over g(Ar, A^) the values of A[u}), eq.(2), of each dressed 
TLS. In amorphous insulators, like amorphous silica, there 
are TLSs with A,, up to about oj* =5 K [24]. The up- 
per cut-off of the distribution is usually larger than the 
frequencies of the flexural modes of interest, lOco < w*. 
Integrating over A^ and Arem we find that the density 
of TLSs per unit volume and unit energy available for 
direct (resonant) excitation processes is given by P, in 




Fig. 3: Sketch of the spectral function of a TLS coupled to a 
sub-ohmic bath. Dashed line shows the off-resonant contribu- 
tion. Dot-dashed line shows the main broadened peak. For 
comparison, the thick vertical line shows the spectral function 
of a non- interacting TLS. 

agreement with the known result that TLSs in amorphous 
systems give rise to a finite density of states at low ener- 
gies [25]. In addition, we find a contribution coming from 
the non-resonant part of the spectral fimction of each TLS, 
^tot~^'^''i^) ~ 2Pab\/uJco/(^- The divergence as w ^ 
arises from the A~^^ dependence of A{iLj), eq.(7). 

Dissipation due to the TLSs. — We assume that 
the externally excited flexural mode of interest, (fco,a;o), 
is linearly coupled, with the coupling constant shown in 
the derivation of the hamiltonian H, 7(Ag/Ao), to a con- 
tinuum of excitations whose spectral strength, Atot{uj), 
is given by the summation over all TLSs of the function 
A{uj) calculated for each one. The ratio (Ag/Ao) can be 
approximated as 1, due to the negligible role played by 
strongly biased TLSs. 

The transition rate of the mode ko occupied by n 
phonons to the mode with n — 1 is calculated from Fermi's 
golden rule, and the energy loss per cycle and unit volume 
AE of the mode will correspond to this transition rate 
multiplied by the energy of a phonon Hloq and the period 
2Tr/uJo: 

27r 27r / fc^ \ ^ 

AE~ — xhuJoX—n[ 7-^ AtoM (8) 

ujQ fi. V V^J 

where Afot(w) = P + Paby/oJco/oJ is the sum of the res- 
onant and non resonant contributions arising from inte- 
grating over the distribution of TLSs, as discussed in the 
preceding paragraph. 

The inverse quality factor Q~^{uJo) is given by 
Q^^{ujo) = AE/2nEo, where Eq is the energy stored in 
the mode per unit volume, Eq ~ nhujo/twL, leading to 
the following expression at zero temperature: 

1/4 

Q-i(wo) ^ lOi^/^^ i-j afaV^Aot(wo) (9) 

Experiments are done at finite temperatures, and, some- 
times, in systems where the oscillator is driven strongly 
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out of equilibrium. We take these effects into account by 
calculating Atot{^,T) in the presence of a finite distribu- 
tion of excited vibrations, which can include a non-thermal 
contribution. If we only keep one phonon processes in the 
calculation of Atot(w, T), as in the previous discussion, the 
subtraction of absorption and emission processes cancel 
the temperature dependence on the number of phonons. 
The only temperature dependence is due to saturation ef- 
fects in the absorbtion properties of the TLSs when their 
environment contains many quanta of energy Aren- The 
final expression for the inverse quality factor at finite tem- 
peratures is: 



E 



1/4 



Ptanh 



Potb 



Wo 



(10) 



Until now prevalence of one-phonon processes in the in- 
teraction among TLSs and vibrational modes has been 
assumed, but at temperatures much higher than the fre- 
quencies of the relevant phonons, multi-phonon processes 
need to be taken into account. We include this effect 
assuming overdamped dynamics for the TLSs, so that 
A(w,A,e„,T) = T(A,e„,T)/(l + [wt(A„„,T)]2), where 
T(Arera,T) = r^^(Aren,T), and making use of the rela- 
tion between Q~^{uj,T) and Aren,T) usually found 
in the context of the standard tunneling model approach 
to disordered bulk systems [24-27,38], which in our case 
translates into: 



cosh 



OJT 



1 + (ujTy 



(11) 



Here e = \/A^g„ -|- A^ and u — Aren/e, with Umin fixed 
by the time needed to obtain a spectrum around the 
resonance frequency of the excited mode and ^max esti- 
mated to be at least of the order of 5 K [24]. The lim- 
its of integration must be such that only the overdamped 
TLSs are included. A given TLS is overdamped when 

A^en < r(A^e„,r) ^ A^e« < [3006^^^] . HcnCC, 

we obtain for the contribution to the inverse quality fac- 
tor from overdamped TLSs: 



(12) 



Therefore, in a range of energies a;/„„ci < ^ ^co, the at- 
tenuation coming from these TLSs shows a Q~^{T) ~ T'^/'^ 
dependence, in qualitative agreement with the experiment 
[18] on Si nanobridges. The total inverse quality factor is 
the sum of eq.(12) plus the temperature independent con- 
tribution arising from off-resonant processes induced by 
underdamped TLSs, eq.(lO). This last equation must in 
principle be corrected by taking into account the decrease 
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4: Fittings of the result of eq.(13) to experimental data 



in the number of underdamped TLSs as the temperature 
is raised, but it is a weak effect and will be neglected, lead- 
ing to the final expression for the total attenuation of a 
flexural mode: 



(13) 



where. 



C(c^o) 



i3/2^£;9/4 



,-1/2 



150P7V^^ r(i + '^)(i-2j^) 



t^w^/^E*uJo 



4/3 



(14) 



As only part of the beam is amorphous, P is to be replaced 
by P ■ Vamorphous/^4otai- To dcscribc the resuhs in [18], 
as shown in fig. [4], we assumed P ■ Vamorphous/Vtotai 
10^'*J~^m~3^ compatible with P values reported in amor- 
phous glasses [24], and 0.1 < Vamorphous/Vtotai < 1- The 
slope of the T^/^ contribution gives the value of 7, used 
as fitting parameter. We obtain 7 ~ 5 — 10 eV, which is a 
reasonable value [39,40]. There are two limitations on the 
range of applicability of our results. The increasing role of 
interactions between TLSs as T is lowered [41], and pos- 
sible cooperative effects when the mode is strongly driven 
[16] cause deviations, which are manifested in the satu- 
ration observed in fig. (4), not explained by our fit, which 
predicts Q^^ ~ 5 • 10"''. On the high-temperature side, 
a point is reached when the rate r~^(Aren,r) changes to 
an Arrhcnius-like behavior [24]. 

Conclusions. — We have studied the damping of me- 
chanical oscillations in nanoscopic devices due to their in- 
teraction with TLSs. This coupling is the main mecha- 
nism of relaxation of phonons in disordered insulators. We 
have analyzed the changes induced in the spectrum and 
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distribution of TLSs due to their interaction with the low 

energy oscillations of nano-mechanical devices. Flexural 
modes, with a high density of states at low energies, lead 
to sub-ohmic damping, which can modify significantly the 
distribution of TLSs. The problem of a TLS interacting 
with a sub-ohmic environment is interesting in its own 
right [28,29,36,42-48], and the systems studied here pro- 
vide a physical realization. 

We obtain a temperature independent contribution to 
the inverse quality factor, of a flexural mode, which 
arises from resonant excitations of TLSs, and off-resonant 
processes involving underdamped TLSs. We find, in addi- 
tion, a contribution which increases as T^/^, arising from 
overdamped TLSs. The off-resonant contributions imply 
that the externally excited vibration loses its energy to 
TLSs which, in turn, decay into other acoustic modes. 
Hence, off-resonant contributions can only be present if 
the number of thermally excited modes is large, a situa- 
tion fulfilled in most present experiments. 

We have made mimerical estimates for the expected dis- 
sipation for a few representative devices. We have assumed 
that a fraction of the device shows amorphous features, 
and contains a distribution of TLSs similar to that found 
in amorphous insulators. The main uncertainties in our 
calculation are due to the lack of information on the TLSs 
distribution, the coupling strength, and the fraction of the 
total volume of the device that they occupy. Decreasing 
volume and number of modes may lead as well to fluctua- 
tions around our predictions, which use continuum distri- 
butions. 
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